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In this paper, which consists of two parts, an extension of the magneto-ionic theory is 
systematically developed by employing the linearized, single fluid, continuum theory of 
plasma dynamics. The cubic equation specifying the square of the phase velocity of a plane 
wave, whose direction of propagation makes an arbitrary angle with that of the static ex- 
ternal magnetic field, is derived. In the first approximation in which only the phase velocities 
of the order of the velocity of the electromagnetic waves in free space or higher are sought, 
the general cubic equation is found to degenerate into a variant of the well-known Appleton- 
Hartree equation. The regions of the normalized source frequency 12 and the normalized 
gyromagnetic frequency R, in which the ordinary and the extraordinary, purely transverse 
electromagnetic modes propagate are discussed with the help of a construction in the il 2 -R 2 
plane and with emphasis on the procedure used in the designation of the two modes. It is 
found that the dispersion curves for the arbitrary direction of propagation do not continu- 
ously go over to those for the direction parallel to that of the static magnetic field. An 
inconsistency in the designation procedure for the modes is found to exist in the magneto- 
ionic approximation. 



1. Introduction 

A vast majority of the investigations on the char- 
acteristics of wave propagation in the ionosphere are 
based on the application of the magneto-ionic the- 
ory. An important assumption of the magneto-ionic 
theory is that the plasma of the ionosphere is incom- 
pressible, with the result its macroscopic properties 
are adequately described by considering it to be a 
dielectric medium characterized by a tensor dielec- 
tric constant. The magneto-ionic theory has been 
successful in explaining the variety of wave phenom- 
enon in the ionosphere, especially in situations where 
the effect of the alternating electric field predomi- 
nates over that of the pressure gradient. This suc- 
cess is to be expected since, in the above circum- 
stances, the dielectric description of the plasma is 
a reasonably good approximation. When the effect 
of the pressure gradient is not negligible, the mag- 
neto-ionic theory fails and consequently requires 
extension. The extension of the magneto-ionic the- 
ory by employing a linearized, single fluid contin- 
uum theory of plasma dynamics is systematically 
developed in this paper by considering the plasma 
of the ionosphere to be an unbounded, collision-free, 
and macroscopically neutral mixture of gas of mobile 
electrons and stationary ions. 

At the outset, a brief discussion on the notation 
used in this paper is in order, especially since it dif- 
fers from the conventional notation employed in the 
ionospheric literature [Ratcliffe, 1959;Budden, 1961]. 
In problems concerning the wave propagation in the 
ionosphere, three frequencies occur; namely, (i) the 



source frequency, co, (ii) the plasma frequency, co,,, 
and (iii) the gyromagnetic frequency, co c . In the 
ionospheric literature, the following normalized 
parameters are in common usage: 



X- 



-^ and F=^. 

or co 



In the analysis of the properties of wave propaga- 
tion in the ionosphere, the determination of the prop- 
agation characteristics explicitly as a function of the 
source frequency, co, is almost always desired. Such 
a determination is rendered unnecessarily cumber- 
some by adhering to the normalized parameters X 
and F, since the source frequency is mixed in both 
of them. Therefore, instead of using the conven- 
tional notation, the following normalized frequencies 
are used: 



_ CO ! T1 w c 

S2=— and R=— 



In this notation, 12 is the normalized source fre- 
quency and E is the normalized magnitude of the 
external magnetic field. The use of the normalized 
parameters, 12 and R, has the additional advantage 
of effecting considerable simplification of the coeffi- 
cients in the dispersion equation, as will be seen in 
the next section. 

There have been previous investigations of the 
nature of plane wave propagation in an unbounded, 
homogeneous, compressible plasma with an external 
magnetic field. These investigations are primarily 
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restricted to two special directions of propagation; 
namely, those along and across that of the static 
magnetic field. An elegant presentation of the 
results for these two special cases may be found in 
the recent book by Allis, Buchsbaum, and Bers 
[1962]. In this paper, first, a brief derivation of the 
dispersion equation for the propagation of a plane 
wave at an arbitrary angle to the direction of the 
static magnetic field is given. The dispersion 
equation, which is essentially a cubic, is found to 
factor out into a linear and a quadratic equations 
for the special cases of propagation along and across 
the static magnetic field. These two special cases 
are briefly reviewed in this paper, partly for the 
sake of completeness and partly to emphasize the 
new construction in the Q 2 —R 2 space, which is used 
in elucidating the regions of 12 and R in which the 
two modes propagate. This new construction serves 
to localize the regions in the 12— R space where the 
validity of the magneto-ionic theory is of suspect. 

It is found that, except in a boundary layer about 
the direction of the static magnetic field, in the first 
approximation wherein the phase velocity of the 
order of the velocity of the electromagnetic waves in 
free space, c or higher, is sought, the exact dispersion 
equation degenerates into a quadratic equation 
which is just a variant of the well-known Appleton- 
Hartree equation. The analysis of the above 
quadratic equation shows that the phase velocity is 
everywhere of the order c or higher, except in the 
vicinities of the frequencies 12 = 0, 12 3 , and 12 4 , where 
it goes to zero. It is shown that 12| and O4 are, 
respectively, the smaller and the larger of the two 
roots of the equation 12 4 — 12 2 (1 +R 2 ) +R 2 1 2 =0, where 
l=cos d and 6 is the angle which the propagation 
direction makes with that of the static magnetic 
field. Therefore, the results of the first approxima- 
tion may not be valid in the neighborhood of the 
three frequencies, 12=0, 12 3 , and 12 4 . 

A closer examination of the exact dispersion 
equation for 12=0 shows that the validity of the 
results of the first approximation is not impaired 
even though the phase velocity goes to zero at that 
frequency. 

In the second approximation wherein the phase 
velocity of the order of the acoustic velocity a, in the 
electron gas or lower, is sought, the exact dispersion 
equation reduces to a linear equation. The study 
of the above linear equation shows that the dispersion 
curve has two branches. The first branch starts 
with a value of infinity for 12=12 3 , reduces progres- 
sively as 12 is increased, and goes to zero for Q—RL 
The second branch starts again with a value of 
infinity for 12=12 4 , reduces progressively as 12 is 
increased, and attains asymptotically the value a 
for 12=oo. Obviously the results of the second 
approximation are also not valid in the neighborhood 
of the frequencies 12=12 3 and 12 = 12 4 . 

A third approximation wherein the phase velocity 
is of the order ^ac~ is, therefore, sought, and the 
exact dispersion equation degenerates then into 
another quadratic equation. The analysis of this 
quadratic equation yields finite, nonzero values for 



the phase velocity for 12=12 3 and 12 4 , and these phase 
velocities are of the order -y/aFo and hence are an 
excellent approximation for the phase velocities at 
12=12 3 and 12 4 . Moreover, it is found that the phase 
velocities obtained in the third approximation merge 
with those given by the first and the second approxi- 
mations, as the frequencies are sufficiently removed 
from the values 12 3 and 12 4 , to the lower or the higher 
sides, respectively. From a synthesis of all the re- 
sults given by the first, the second, and the third 
approximations, emerge the complete dispersion 
curves which are valid for all directions of propaga- 
tion except very close to that of the static magnetic 
field. In all, there are three modes and these have 
been named as the modified plasma mode, the modi- 
fied ordinary electromagnetic mode, and the modified 
extraordinary electromagnetic mode, respectively. 
Only the modified extraordinary electromagnetic 
mode has two branches. 

In the axial boundary layer, which is a narrow 
region in which the propagation direction is very 
close to that of the static magnetic field, the dis- 
persion curves undergo very rapid changes as a 
result of the decoupling of the plasma mode. Con- 
sequently, dispersion curves which are obtained for 
the arbitrary direction of propagation to that of the 
magnetic field do not uniformly go over to those for 
the propagation direction along that of the static 
magnetic field. It is found that the dispersion in 
the axial boundary layer undergoes rapid changes 
in the neighborhood of the two frequencies 12 =R 
and 12=1. By using boundary layer techniques, the 
dispersion equations in the boundary layer in the 
close neighborhood of the two frequencies 12 =i? and 
12=1 are derived and analyzed. This analysis clari- 
fies the manner of coupling between the longitudinal 
plasma waves and the transverse electromagnetic 
waves for the case of propagation in the directions 
other than that of the static magnetic field. When 
the propagation direction is changed from that of 
the magnetic field, for the case B 2 <1, the ordinary 
electromagnetic mode couples into the plasma mode 
and vice versa, in the close neighborhood of the 
frequency 12=1, and when the phase velocity is 
approximately equal to c . For the case 2? 2 >1, in 
addition to the above coupling, the plasma mode 
couples into the extraordinary electromagnetic mode 
(Whistler mode) and vice versa, in the close neighbor- 
hood of the frequency 12=1 and when the phase 
velocity is slightly lower than c . Also in the neigh- 
borhood of the frequency 12= R, when both the extraor- 
dinary electromagnetic mode and the plasma mode 
are slowed down to approximately the acoustic ve- 
locity a in the electron gas, they couple back into 
each other. The smaller and the larger of the two 
frequencies 12 = 1 and 12 = 7? are found to be respec- 
tively the approximations to the frequencies 12=12 3 
and 12=12 4 in the axial boundary layer. As the 
propagation direction is changed so as to approach 
that perpendicular to the static magnetic field, the 
only essential changes in the dispersion are the de- 
crease of 12 3 and the increase of 12 4 . 

The general dispersion equation is briefly derived 
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in section 2. In sections 3 and 4, the two special 
cases corresponding to the direction of propagation 
along and across that of the static magnetic field are 
discussed. The first order approximation to the 
exact dispersion equation is treated in section 5. 
The higher order approximations for the region ex- 
ternal to the axial boundary layer as well as the 
dispersion in the boundary layer are contained in 
part II of this paper. 

2. Formulation of the Problem 

Consider an unbounded, homogeneous plasma 
which, for the sake of simplicity, is assumed to be a 
lossless and macroscopically neutral mixture of gas 
of mobile electrons and immobile ions. It is proposed 
to restrict attention only to the linear, time-harmonic 
problem; the harmonic time dependence of the form 
e -i<at j s implied f or a u the field components. It is 
further assumed that the drift velocity of the electrons 
is zero so that the plasma as a whole may be consid- 
ered as stationary. A uniform magnetic field B is 
assumed to be impressed externally throughout the 
plasma region in the ^-direction, where x, y, and z 
form a right-handed rectangular coordinate system 

(fig.i). 

Let N be the average number density, P the 

-> 
pressure deviation from the mean, and V the velocity 

-> -» 

of the electrons. Let E and H be the alternating 
electric and magnetic fields. The linearized time- 
harmonic hydro dynamic equation of motion for the 
electrons is 



-iumN V=N e(E+VXzB )-vP 



(1) 



where e is the charge and m is the mass of an election . 
The equation of continuity after being linearized and 
combined with the equation of state is given by 



a 2 mN V'V=io)P 



(2) 



where a is the velocity of sound in the electron gas. 
In addition, the electric and the magnetic fields 
satisfy the following time-harmonic Maxwell's equa- 
tions 

VXS-icW? (3) 



\vXH=-ia>e E+N eV 



(4) 



where ju and e are the permeability and the dielec- 
tric constant of free space. 

It is desired to investigate the characteristics of a 
plane wave propagating in the plasma medium. 
-> , A -> 

Let k be the propagation vector such that zXk 
coincides with the ?/-axis and let 6 be the angle that 
-> 

k makes with the direction of the static magnetic 
field, so that k x =k sin d=kn, k y =0 and k z =k cos 
d=kl. Therefore all the field components will have 



Figure 1. Geometry of the problem. 



the spatial dependence of the form e ik(nx+lz) so that 



■^— =ikn. ^r~=0; 

dx dy 



dz 



= ikl. 



(5) 



The elimination of P from (1) with the help of (2) 
and the use of (5) leads, after some rearrangement, to 
the following simultaneous equations specifying 
V x , V y , and V z in terms of E x , E V) and E z \ 

( X J^l\ Vi _kWnl y Je_ E _p_ E 
\ ora / ara coma arma 



ico r k 2 a 2 n 2 



v x +v v 



ico c k 2 a 2 nl 



V z =- 



eo) c 



K 



ze 
coma 



E„ 



¥a 2 nl 



V x +(l-^)v^E lt (6a,b,c) 



where 



and 



a=l-- 



(7) 



eBo 



--—■ — - is the gyromagnetic frequency of the 

electrons. The substitution of V x , V y , and V 2 
obtained from the solution of (6) into (4) enables 
the right-hand side of (4) to be rewritten as 

—icoe^e-E, where e" is the dielectric tensor. Then, 
from (3) and (4), it is found that 



vxvxe= 



772 € ^ 

^0 



(8) 
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where C : 



1 



Moeo 



is the velocity of electromagnetic 



waves in free space. On writing (8) in component 
form, the following set of equations specifying E x , 
E y , and E z is obtained: 



-iD 2l 

. Ai 

where 



iD l2 
D 22 

_k 2 l 2 C 2 



A." 




-E- 


-iD 23 




Ey 


A 3 _ 




-E 2 _ 



--[D][E}=0 



(9) 



Vn -^r ~ 1+ mdv — it) 

a, 2 , w c A £W\ 
oj 2 coaA\ a) 2 / 



A 2 =A,=- 



D»=D„ = - 



FlnCl , co 2 Fa 2 »Z 



VD 



n _F(7g , , a, 2 / &V\ 



n — n — v 

JJ 2Z — Urn — 7 



co c k 2 a 2 nl 
aaD co 2 






(10) 



and o) p is~the plasma frequency of the electron and is 
given by ca 2 p = — — . The determinant of D should be 

TYl €q 

equal to zero in order that (9) may have a nontrivial 
solution and this condition, together with the help of 
(10), after considerable manipulations leads to the 
following cubic equations in X 2 : 



\ Q + A 1 \* + A 2 \ 2 + A 3 = 



(11) 



where X=ir is the phase velocity of the plane waves 

propagating in the medium. The expressions for 
the coefficients A u A 2 , and A z in (11) may be simpli- 
fied to yield the following: 



A=- 



U 2 Cl 



A 



l{2&-2Q 2 (2+R 2 )-\-2+R 2 (l + l 2 )} 
+8{2 i -a 2 (2+l 2 R 2 )+l})=^ 



Ai= &C% [ 4_fi 2 (l+^ 2 )+i? 2 Z 2 



A 



+2n 2 8{n 2 -i-i 2 R 2 }]=4^ 



(12a) 



(12b) 



A* 



&Cta 2 (W-l 2 R 2 )_A,, 



A =(a 2 - 

where 



A A 

-l)[Q 4 -0 2 (2+£ 2 ) + l] 



x a o 



-±,B= 






(12c) 
(12d) 

(13) 



3. Propagation Along the Static Magnetic 
Field 

Before proceeding to obtain the characteristics of 
plane waves propagating at an arbitrary angle 6 with 
respect to the direction of the external magnetic field, 
the characteristics of plane waves propagating along 
and across the static magnetic field will be briefly 
reviewed. For propagation along the magnetic 
field, 0=0 and the dispersion equation (11) may be 
factored into the following two equations: 



and 



where 



Br- 



x 2 - a2Q2 -0 



\ i +B 1 \ 2 +B 2 =0 

2a 2 C 2 (Q 2 -l-R 2 )_ 



B 2 = 



B 
Q 2 Ct(&-R 2 )_ 



B w 
~B 



(14a) 



(14b) 



(15a) 



Bo, 



(15b) 

(15c) 



-B B 

.B =[Q 4 -fl 2 (2+£ 2 )+l]. 
The phase velocity obtained from (14a) is given by 



v -=w- 



•1 



(16) 



The mode whose phase velocity is given by (16) 
propagates only for £2>1 and has a phase velocity 
equal to the acoustical velocity in the electron gas in 
the limit of infinite frequency. This mode is 
obviously the plasma mode. Note that (16) is 
independent of the static magnetic field and therefore 
the same as in an isotropic plasma. This is to be 
expected since the static magnetic field does not 
exert a force in a direction parallel to itself. 

The dispersion equation (14b) is independent of a 
and hence is the same as in an incompressible plasma. 
It may be solved using a procedure to be described 
in connection with the case of propagation at an 
arbitrary angle and the resulting phase velocities of 
the ordinary and the extraordinary electromagnetic 
modes are given by 
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In figure 2a, the regions of propagation of the ordi- 
nary and the extraordinary modes are indicated. 
The ordinary electromagnetic mode is found to prop- 
agate for £2>Oi and the extraordinary electromag- 
netic mode for £2<7? and 12>12 2 , where 



12 1( 2 



*i+\- 



/(!)' 



+1. 



(18) 



An examination of figure 2a reveals a fact which 
does not appear to have been emphasized previously 
in the literature. For i? 2 >l/2, there is at least one 
mode propagating in the entire frequency range. 
But for R 2 <^l/2, there is a frequency band which lies 
entirely below 2 = 1 and for which neither mode 
propagates. As R 2 is reduced from the value 1/2, 
this cutoff band increases progressively until finally 
for the isotropic case, this cutoff band extends from 
tt=0 to 0=1. 

The phase velocities of the three modes, given 
respectively in (16) and (17a, b) are sketched in 
figure 2b for the case /^ 2 <1 and in figure 2c for the 
case i? 2 >l. The ordinary electromagnetic mode is 
seen to have a phase velocity which starts with a 
value of infinity for 12= 12 1? becomes smaller as Q is 
increased and attains asymptotically the free space 
value Co for infinite frequencies. The higher fre- 
quency branch of the extraordinary mode has a 
phase velocity which starts with infinity for 12=12 2 
and reduces asymptotically to the value C for infinite 
frequencies. The low frequency branch of the 
extraordinary mode, usually referred to as the 
"whistler mode/' has a phase velocity which starts 




Figure 2a. Regions of propagation of the ordinary and the 
extraordinary EM modes for propagation along the static 
magnetic field. 




Figure 2b. Dispersion curves for propagation along the 
direction of the static magnetic field '(R 2 = }£). 

with a value zero for £1=0, increases and reaches a 
maximum value for &=R/2 and for further increase 
in 0, the phase velocity reduces and finally becomes 
zero for tt = R. The maximum value of the phase 
velocity which is attained for tt=R/2 is equal to 
C R(R 2 +4:)~ l/2 which is always less than C but 
approaches that value as R is increased. Only in 
the frequency range 0<Q<R/2, the phase velocity 
increases with frequency. In figure 2a, the region to 
the left of the critical line R 2 =m 2 , is the region of 
the "whistler mode," where the phase velocity in- 
creases with the frequency. 
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Figuke 2c. Dispersion curves for propagation along the 
direction of the static magnetic field (Ji 2 = ji). 

4. Propagation Across the Static Magnetic 
Field 

For propagation across the magnetic field, 6=11/2 
and the general dispersion equation (11) may be 
factored into the following two equations: 



where 



w 



o 



and 



2 -l 
K+B 1 \ 2 +B 2 =0 



(19a) 



(19b) 



B 1 - 



B 2 = 



B 10 _ [V 2 C 2 (n 2 - l-R 2 ) + Q 2 a 2 (tf-l)] 
B 

B 20 a 2 C 2 W 



B, 



(20a) 



(20b) 



Bo Bo 

Bo=W-tt 2 (2+R 2 )+l. (20c) 

The phase velocity obtained from (19a) is given by 

C n 



Vo 



V^ 2 -i 



(21) 



The mode whose phase velocity is given by (21) prop- 
agates only for Q>1 and has its phase velocity 
asymptotically approach the free space electro- 
magnetic wave velocity. Also this mode is the same 
as in an isotropic plasma. Evidently this mode is 
the ordinary electromagnetic mode with its electric 
vector parallel to the direction of the static magnetic 
field. 

The dispersion equation (19b) contains both a 
and Co and hence pertains to modes obtained as a 
result of mixing of the extraordinary electromagnetic 
wave and the plasma wave. These modes will be 
designated as the modified extraordinary electro- 
magnetic mode (MEM[x]) and the modified plasma 
mode (MP) and will be denoted by the subscripts 
mx and mp respectively. On solving (19b) and 
simplifying the resulting expressions by exploiting 
the fact that a 2 /Cp<<l, the following expressions 
for the phase velocities of the modified extraordinary 
electromagnetic mode and the modified plasma mode 
are easily obtained: 

VC 2 M 2 (Sl 2 — S]i)T /2 „ 
I^=[ ° jg 3) \ forG 2 <Q<«> (22a) 



where 



Vt=l- 



R 2 



a+i? 2 



1+ 



ci 



It may be easily established that fi?<&3<122. 
V, 1 ,=[ ^ Q2( ^~ 0g) ] ' /8 for S2 1 <0<fi 3 



(22b) 

Also, 
(23a) 



= aV^[§ (£2 2 -«l)+(-^o) 1/2 J 

for|Q-ft 3 |<0(10- 4 ) (23b) 



V m 



r d 2 o 2 "| 1/2 

1(8 2 -GDJ 



for 12 3 <^fi<^oo. 



(23c) 



The phase velocity of the modified extraordinary 
electromagnetic mode starts with the value infinity 
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for Q=ti 2 , becomes smaller as 12 is increased and 
asymptotically approaches the free space electro- 
magnetic wave velocity C in the limit of infinite 
frequency. The phase velocity of the modified 
plasma mode starts with the value of infinity for 
Q=ti u becomes smaller as il is increased and reaches 
the value C () approximately for 12=1. At about 
12=12 3? the phase velocity rapidly falls and reaches 
the value of the acoustic velocity in the electron gas 
in the limit of infinite frequency. It is not difficult 
to verify that the modified plasma mode is a pre- 
dominantly transverse wave for 12!<^<^ 3 and a 
predominantly longitudinal wave for 12 3 <12<oo. 
Obviously the coupling between the transverse and 
the longitudinal waves takes place in the close 
neighborhood of 12=12 3 . 

In figure 3a, the regions of propagation of the 
modified extraordinary electromagnetic mode and the 
modified plasma mode are indicated. Also the phase 
velocities given by (21), (22a), and (23) are sketched 
in figure 3b for the case II 2 =1/2. 



MODIFIED PLASMA MODE 




Figuee 3a. Regions of propagation of the modified plasma and 
the modified extraordinary EM modes for propagation across 
the static magnetic field. 



5. Propagation at an Arbitrary Angle to the 
Direction of the Static Magnetic Field 
(First Approximation) 

Except for the cases of propagation along and 
across the static magnetic field, the dispersion equa- 
tion (11) which is a cubic in X 2 , does not exactly factor 
out into two simpler equations. Nevertheless, it is 
possible to exploit the fact that a 2 /Co<<l and 
perform a perturbation analysis of (11), which in 
essence leads to a separation of (11) into simpler 
equations which are either linear or quadratic in X 2 . 

The free space electromagnetic wave velocity C 
is approximately equal to 3X10 8 m/sec and the 




Figure 3b. Dispersion curves for propagation across the 
direction of the static magnetic field (R 2 = }{>). 



acoustic velocity in the electron gas is of the order 10 4 
m/sec and may be taken for convenience to be equal 
to 3X10 4 m/sec. Hence 5 = a 2 /C 2 =W- s which is 
negligible in comparison with unity. It is conven- 
ient to rewrite (11) together with (12) as follows: 



<r 



6 +\4 [A (1) +£s ^ <2) ]+x 2 \_AP+^i < 



'•] 



M.=0 

(24a) 
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where 



A? 



n 2 O 2 



A 



[2Q4-2fi 2 (2 + # 2 )+2+£ 2 (l + / 2 )] = 



2 f7 4 A <« 

=^P [fi 4 -8 2 (14-i? 2 )+iJ*Z*]=^2. 



=A=- 



4 Cfe 2 (0 2 -Z 2 i? 2 )_^> 



A 



A 



A) 



and 



A 



.a 2 c* 



A 



[2fi 2 (G 2 -l-Z 2 .R 2 )]. 



id) 
_£iio_ 

~ A 

(24b) 

(24c) 

(24^ 
(24e) 

(24f) 



If the phase velocity is of the order C 0) it is seen from 
(24a) that there are three terms of the order Co and 
the other three are of the order C%(a 2 /Cl). Hence 
only the terms of the order Co in (24c) need be re- 
tained resulting in the following quadratic equation 
in X 2 : 



\*+A™\*+A?>-. 



= 0. 



(25) 



It is to be noted that on setting a equal to zero in 
(11) and (12), (25) results. This is equivalent to 
saying that had the magneto-ionic theory been used 
from the outset, the dispersion equation (25) will be 
obtained. It is therefore obvious that (25) is a 
variant of the well-known Appleton-Hartree equa- 
tion. It is proposed to analyze the dispersion equa- 
tion in a more systematic way than has hitherto been 
done and then proceed to determine the perturbation 
caused by taking into account the compressibility of 
the medium. 




Figure 4. Regions of propagation of the ordinary and the 
extraordinary EM modes for the arbitrary direction of propaga- 
tion (first approximation) . 



The solutions of (25) are given by 



A (1) 1 l/A a) \ 2 

^-k^hli 4 ?) - A - A - (26) 

The lower and the upper signs in (26) correspond 
respectively to the ordinary (0) and the extraor- 
dinary (e) electromagnetic modes. When 



I 1 l/A {1) \ 2 



> 



Ai 



2A 



and hence the solution (26) which corresponds to the 
positive sign in front of the radical will yield a posi- 
tive value for X 2 resulting in the propagation of the 
corresponding mode. In order to find out the region 
in the tt 2 —R 2 space in which A { 2oA <C0, the critical 
lines corresponding to ^4o=0 and A2o=0 are drawn 
in figure 4. The critical lines corresponding to 
A =0 are given by tt 2 =l and # 2 =iT 2 (ft 2 — l) 2 and 
that corresponding to A { io are given by fi 2 =0 and 
R 2 =Sl 2 (ti 2 -l)/(ti 2 -l 2 ). It may be easily shown 
that the line R 2 =tt 2 (tt 2 —l)/(2 2 -l 2 ) is above the line 
B 2 =Q- 2 (Q 2 -l) 2 ior fi>l. The signs of A and 
A ( io in the various regions into which the ti 2 —R 2 
space is divided by the above mentioned critical 
lines are indicated in figure 4. In the regions 
marked 1, 2, and 3, A ( ^A <^0. In the regions 1 
and 3, A <^0 and hence, the propagating mode will 
correspond to the lower sign in (26) and is the 
ordinary electromagnetic mode. This is indicated 
by shading with the horizontal lines. In the 
region 2, A o ^>0 and the propagating mode will 
correspond to the upper sign in (26) and is the 
extraordinary electromagnetic mode which is indi- 
cated by shading with the vertical lines. 

Using (24a, b) and (12d), it can be shown that 



v 4f ) -ASA. 



--C*Q 2 R 2 



(ft 2 -l) 2 / 2 



Q 2 T? 2 "I 

-^f-(l-H 2 J>0 (27) 



and hence X 2 ,,*, are always real. If A 2 l o ^4 >0, 
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and hence the signs of both X 2 , and X 2 , will be the 
same as A ( $/A Q . It is therefore evident that both 
the ordinary and the extraordinary electromagnetic 
modes will propagate if A ( $/A Q <iO. I n figure 4, 
the critical lines corresponding to A =0 and A ( 2 o = 
0, and also the critical lines 12 2 = and R 2 =(tt 2 — 

l)/( 12 2 — — - — ) corresponding to A ( io=0 are drawn. 

For R 2 yo the line R 2 =Q 2 (V 2 -l) 2 /(n 2 -^^\ may 

be shown to be always above the line R 2 =tt~ 2 
(ft 2 -l) 2 . Also the line J? 2 = Q 2 (fi 2 -l)/(12 2 -/ 2 ) 
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found to be above the line 7i >2 =(S] 2 -l) 2 /l2 2 -^t^") 



2 

In the 



for S2 2 >1. Furthermore, l 2 <^j^<l. 

regions for which A$A ^>0, the signs of both 
A ( io and A are indicated in figure 4. It is found 
that A { ii/A <C0 in the regions marked 4, 5, and 6 
and, therefore, both modes propagate in those 
regions. The line R 2 = constant intersects the line 
# 2 =£} 2 (ft 2 -l)/(ft 2 -/ 2 ) for two values of SI 2 given by 



l+R 2 



=F 



ve 



+R 2 \ 2 
2 )~ 



■R 2 l 2 . 



(28) 



In (28), 12§ a nd Ql correspond to the upper and the 
lower signs respectively and are always real. Also 
the line R 2 = constant intersects the line R 2 =ti~ 2 
(£2 2 — l) 2 for two values of ft 2 , namely ft? and 0% 
whose expressions are given in (18). From figure 4, 
it is clear that the ordinary electromagnetic mode 
propagates in the two frequency ranges 0<0 2 <12i 
and 1<12 2 < °°^and the extraordinary electromagnetic 
mode in the ranges ^ 2 <fi 2 <04 and 122<12 2 <co. 
For the sake of convenience in the redesign a lion of 
these modes after the modification resulting from the 
inclusion of the effect of the compressibility of the 
medium, the ordinary electromagnetic modes propa- 
gating in the frequency ranges 1<12 2 <co and 
0<^ 2 <C^3 are denoted by 0\ and U , respectively, 
and the extraordinary electromagnetic modes propa- 
gating in the frequency ranges 122<£2 2 < °° and 
£2i<0 2 <^4 are denoted by e 1 and e lI 9 respectively. 

With the help of figure 4, at least one mode is seen 
to propagate in any frequency range, for R 2 ^>1/1 2 
(l+l 2 ). When IP<l/P(l+F), there is a frequency 
band which lies entirely below 12 2 = 1 and in which 
there is no propagation. As R 2 is decreased from the 
value l/l 2 (l-\-l 2 ), this cutoff band increases progres- 
sively in width until finally for the isotropic plasma 
this cutoff band extends from £2 2 =0 to tt 2 =l. 

As the direction of propagation of the plane wave 
is changed so as to approach the direction normal to 
the static magnetic field, I 2 becomes smaller and as a 
consequence Q% becomes smaller and Qt becomes 
larger. In the limiting case of the propagation 
across the static magnetic field, l 2 =0 and therefore, 
Oj|=0 and tit=l+R 2 . Note that Q? and $ do not 
depend on I 2 and therefore, on the direction of propa- 
gation of the plane waves. Therefore, as the propa- 
gation direction of the plane waves approaches the 
direction perpendicular to the static magnetic field, 
the frequency range of propagation of the ordinary 
electromagnetic mode n decreases and that of the 
extraordinary electromagnetic mode e 11 increases pro- 
gressively until for the limiting case of propagation 
across the static magnetic field, the 11 mode dis- 
appears completely and the propagation range of the 
e 11 mode extends from $<£2 2 <l+# 2 . It is to be 
noted that the frequency bands of propagation of 
both the l mode and the e l mode remain unaltered 
as the propagation direction is changed. 




Figure 5. Variation of the critical frequencies tt^ and S2 4 with I. 



As the direction of propagation of the plane wave is 
changed so as to approach the direction of the static 
magnetic field, I 2 becomes progressively larger and 
reaches the limiting value of unity. As I 2 becomes 
larger, £2r< increases for all values of R 2 in such a way 
as to approach t he value which is equal to t he smaller 
of the two quantities, R 2 and 1, as shown in figure 5. 
Also, 0% decreases continuously for all values of R 2 
in such a way as to approach the value which is 
equal to the larger of the two quantities, R 2 and 1. 
In the limiting case of propagation along the static 
magnetic field, the critical line corresponding to 
Hi for 7» >2 < 1 and that corresponding to Q% for R 2 > 1 
join to form the critical line Q?=R 2 whereas the 
critical line 12 2 = 1 corresponding to that of 12 2 for 
R 2 < 1 and of Q% for R 2 > 1 disappears. It is found 
with the help of (12d) and (26) that when the critical 
line 12 2 =1 disappears, the sign in front of the radical 
in (26) will be reversed for ft 2 <0, with the result in 
the limiting case of propagation along the static 
magnetic field, the U and the e 11 mode will become 
the extraordinary and the ordinary electromagnetic 
modes respectively. Finally the situation depicted 
in figure 2a is obtained and the ordinary electro- 
magnetic mode propagates in the range O^O^ °° 
and the extraordinary electromagnetic mode in the 
two frequency ranges 0<12</t and r 2 2 <12< 00. The 
reason for the discontinuous disappearance of the 
critical line 12 2 =1 for the case of propagation along 
the direction of the static magnetic field is clarified 
when the compressibility of the plasma is taken into 
account and is essentially due to the decoupling of 
the plasma mode, which remains coupled to the 
electromagnetic mode for all directions of propagation 
other than that of the static magnetic field. 

The expressions for the phase velocities of the 
ordinary and the extraordinary electromagnetic 
modes may be written down explicitly with the help 
of (12d), (24b, c), (26), and (27) as follows: 
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(29) 



The upper sign in (29) corresponds to propagating 
extraordinary modes in the frequency ranges 
ftj^ft^flU and 2 <^< oo and the lower sign in 




f\ 5 ,.,44 
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Figure 6a. Dispersion curves for the arbitrary direction of 
propagation in the first two approximations (R 2 =}4). 



(29), to the ordinary electromagnetic mode in the 
frequency ranges 0<£2<£2 3 and 1<12<oo. The 
phase velocities of the two modes as obtained from 
(29) are plotted (solid line) in figure 6a and b for 
for l 2 = 1/2 and for two values of R 2 , namely (i) R 2 =l/2 
and (ii) iZ 2 =4/3. The phase velocities of both the 
1 and the e 1 modes are seen to start with a value of 
infinity for 12=1 and &=tt 2 , respectively, decrease 
continuously as 12 is increased and asymptotically 
approach the velocity C and hence are always 
higher than C . The phase velocities given by (29) 
are obtained from (25) which is an approximation 




0.577 \ 1.414 2 

0.816 



Figure 6b. Dispersion curves for the arbitrary direction of 
propagation in the first two approximations (R 2 = %). 
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to the exact dispersion equation (11) for the case in 
which the phase velocity is of the order C or higher. 
Consequently the phase velocities of the 1 and the 
e 1 modes are excellent approximations to the actual 
values obtained from the solution of (I 1). 

The phase velocity of the n mode stalls with a 
value zero for 12=0, increases as 12 is increased till a 
maximum value of the order of C is reached and 
thereafter decreases for further increase in 12 and 
reaches the value zero for 12^12 3 . Note that in the 
first part of the frequency range 0<C^<^3, the phase 
velocity of the O 11 mode increases with frequency. 
The phase velocity of the e 11 mode starts with a value 
of infinity for 0=12! (which is always less than 1), 
decreases as 12 is increased, attains the value C 
approximately for 12=1 and rapidly decreases to the 
value of zero at £2=12 4 . Since the phase velocity of 
the O 11 mode in the vicinity of 12=0 and 0=12 3 and 
that of e n mode in the vicinity of 12 = 12 4 are con- 
siderably below Cq and since these values are derived 
from (25) which is a valid approximation to t he exact 
dispersion equation (11) only if the phase velocity is 
of the order C or higher, it follows that the phase 
velocity obtained from (25) is of doubtful validity in 
the vicinity of the frequencies 12=0, 12 3 , and 12, and 
therefore a second approximation to (11) is needed 
to examine and obtain the phase velocities in t be 
neighborhood of 12=0, 12 3 , and 12.L This is carried 
out in part II of this paper. 

If a uniform procedure is used in designating the 
modes, it is found that due to the discontinuous 



disappearance of the critical line £2 2 =1 for the case of 
propagation in the direction of the static magnetic 
held, the "whistler mode" has its nomenclature 
changed from extraordinary to ordinary as soon as 
the 1 direction of propagation differs from that of the 
static magnetic field. This inconsistency appears 
only in the first approximation and therefore charac- 
teristic only of the magneto -ionic theory. With the 
higher order approximations to the exact dispersion 
relation and with the effect of the compressibility 
of the plasma thus taken into account, the three in- 
dependent modes which propagate in an electron 
plasma are capable of being designated using a pro- 
cedure which is uniformly valid for all directions of 
propagation. The discontinuous disappearance of 
the critical line 12 2 =1 is shown in part II to be the 
consequence of the decoupling of the purely longi- 
tudinal plasma wave for propagation in the direction 
of the static magnetic field. 
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